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EXISTENCE AND MAXIMAL LP-REGULARITY OF SOLUTIONS FOR THE 
POROUS MEDIUM EQUATION ON MANIFOLDS WITH CONICAL 

SINGULARITIES 

NIKOLAOS ROIDOS AND ELMAR SCHROHE 


Abstract. We consider the porous medium equation on manifolds with conical singularities and 
show existence, uniqueness and maximal L^-regularity of a short time solution. In particular, 
we obtain information on the short time asymptotics of the solution near the conical point. 
Our method is based on bounded imaginary powers results for cone differential operators on 
Mellin-Sobolev spaces and i?-sectoriality perturbation techniques. 


1. Introduction 

The porous medium equation is the parabolic diffusion equation 

(1.1) u'(t) + A(u™(t)) = f{u,t), te(0,To], 

(1.2) u(0) = Mo¬ 

lt describes the flow of a gas in a porous medium, where u is the density distribution. We assume 
here that m > 0 and / = f{X,t) is a holomorphic function of A on a neighborhood of Ran(Mo) 
with values in Lipschitz functions in t on [0,To]. The porous medium equation can be regarded 
as a nonlinear version of the heat equation, which arises for m = 1. In general it can model heat 
transfer, fluid flow or diffusion. While traditionally treated on domains in R", the porous medium 
equation has also been considered on closed compact and noncompact manifolds; we refer for this 
to [3], [U, [IS], [21] and [33]. 

In this article, we study the equation on a manifold with conical singularities. The principal 
application we have in mind is that of a two-dimensional surface with conical singularities, on 
which we imagine the medium as a film. However, both our methods and results extend to higher 
dimensions. We model the manifold with conical singularities by an (n -|-1 (-dimensional compact 
manifold B with boundary, n > 1, endowed with a degenerate Riemannian metric g, which, in a 
collar neighborhood [0,1) x 9B of the boundary 9B, admits the warped product structure 

g{x,y) = dx'^ + x'^h{y). 

Here /i is a (non-degenerate) Riemannian metric on 9B, and {x, y) € [0,1) x 9B. By A we denote 
the Laplacian associated with this metric. It naturally acts on scales of Mellin-Sobolev spaces 
Hp’'^(B), which will be introduced in Section EU At this point it suffices to say that s € R gives 
the local Sobolev regularity with respect to an base space, 1 < p < oo, and 7 € R is a weight 
at the tip. 

We will prove existence and maximal regularity of a short time solution to the porous medium 
equation. Our main tool will be maximal regularity theory for quasilinear parabolic problems, 
based on the properties of the associated linearized problem and a theorem of Clement and Li. We 
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start by considering the Laplacian as an unbounded operator and choosing an appropriate closed 
extension. The possible domains consist of a Mellin-Sobolev space plus a finite dimensional space 
of functions with specific asymptotics near the tip, the so-called asymptotics spaces] see Section 3 
for details. The maximal regularity approach then allows us to examine the effect of the singularity 
on the solution of the problem. 

Finding a maximal regularity extension of the Laplacian is not a trivial task. Certain examples 
of such extensions on Mellin-Sobolev spaces of non-negative order (i.e. s > 0) are known 

explicitly from [53], jUj, |5B]. These will be the starting point of our analysis. We will, however, 
have to go beyond the methods used in these articles: As a first result, we show that there exist 
extensions of the Laplacian of maximal regularity also in Mellin-Sobolev spaces of negative 
order. In fact, for each s € M and suitable 7 , depending on dim(B), we may choose T’(Ag) = 
0 C as a domain in H®’'’'(B) in which the Laplacian has maximal regularity. Here C 
stands for the constant functions. 

The second important step, in order to apply the theorem of Clement and Li, is to control 
the ‘trace space’, i.e. the interpolation space specifying the admissible initial values. This is 
connected to the problem of understanding the domains of the complex powers for cone differential 
operators. For a conic manifold, results on the interpolation between Mellin-Sobolev spaces and 
Mellin-Sobolev spaces plus asymptotics spaces were not known. We manage to obtain suitable 
embeddings for these spaces in our situation. 

The main difficulty then is to show maximal regularity for the linearized term of the problem, 
which is a conically degenerate differential operator of second order. The idea is to first establish 
uniform i?-sectoriality for the operators with frozen coefficients and then to apply perturbation 
arguments in order to construct i?-sectorial local approximations of the resolvent. This allows us 
to show the existence and i?-sectoriality of the resolvent via a Neumann series expression. Our 
first result is the following. 

Theorem 1.1. Let Ai he the largest nonzero eigenvalue of the boundary Laplacian Ag, induced by 
the metric h on 5B. Write 

dimt”) - - A, >0. 

Let 7 satisfy 

dim(B) - 4 f dim(B) - 4 dim(B) 

2 <7<nRn| 

and choose p, q so large that 

dim(B) 2 ^ , dim(B)-4 4 

p q 2 q 

We consider the extension 

(1.5) A, : V{A,) = 7^;+2'T'+2(B) © C ^ 7^;-t'(B). 

Then, for any s > — 1 + + | and for any strictly positive initial value 

uo G (h;+ 2 ’^+ 2 (b) 0 C,H;’^(B))i_, 

there exists a T > 0 such that the porous medium equation (ED, (EDa considered in the space 
L'^{0,T;'Hp''^{M)), has a unique solution 

u G L‘?(0, T; H;+2’'>'+2(B) 0 C) n lTi’«(0, T; '^^’'^(B)). 


(1.3) 


£ = — 


dim(IB) — 2 
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According to Theorem IIL4.10.2 in [5] the solution then automatically belongs to 


^([0, T ]-© C, 


Moreover, we have the following embedding, see Lemma [521 For each e > 0 


n 


s+2-^+e,j+2-j+e, 


iC -A (H 


s+2,'y+2 ^ 




s+2-^-e.7+2-i-E, 


©C. 


The assumptions imply that the domain in (11.51) is a subset of C'®', where tr has to be slightly larger 
than 1. There is not much point here in striving for lower regularity as the specific character of 
the conical singularity would be lost. 

Finally, we analyze the asymptotics of the solution near the tip of the cone as t 0+. To this 
end, we will work with a different extension, namely 


(1.6) A:V(A):=V(A^,)^V(AJ, 

for Ag as above. The point is that we know (D(Ag) rather well in terms of asymptotics and 
regularity. In fact, the asymptotics terms can be computed explicitly from the spectrum of the 
boundary Laplacian Aq. In order to treat the nonlinearity we will have to make the assumption 
that the first nonzero eigenvalue of Aq is < —9/4 for dim(]B) = 2 and < I — dim(B) for dim(B) > 4. 
This implies that the parameter e defined in Theorem II.II is >3/2 for two-dimensional B and > 1 
when the dimension is larger than three. We let p, q, s as before and assume that 


J-1/2 < 7 < min{£ - 2 - ^,0}, dim(B) = 2 

(1-7) j dim(B)-4 ^ ^ ^ jnin { dim(B)-4 ^ g _ ^ dim(B)-2 | ^ dim(B) > 4 ’ 

As a consequence, the first non-constant asymptotics term in ^(A^) will be o{x) near the tip x = 0 
(even o{x^^^) in the two-dimensional case). It is this fact that eventually will allow us to treat the 
nonlinearity in the equation. This approach breaks down in dimension three. 

We obtain the following theorem. 

Theorem 1.2. Let dim(B) ^ 3 and assume that the first non-zero eigenvalue Ai of Aq satisfies 
Ai < —9/4 for dim(B) = 2 and Ai < 1 — dim(B) for dim(B) > 4. Choose s,p, q as in Theorem M.W 
and 7 as in Equation (113. Let A he the closed extension of the Laplacian defined in (11.61) . Then, 
for each strictly positive initial value 

uo€(P(A 2 ),P(AJ)i^,, 

there exists a T > 0 such that the porous medium equation dni). (II3 has a unique solution 

u G L«(0,T;P(A2)) n WF9(0,T;P(AJ). 

This solution automatically belongs to C([0,T]; (T>(A^), (D(Ag)) i : The precise description of 

g 

(D(A^) in Section iTTl below, then gives a clear picture of the possible asymptotics. 

We will start in Section 2 by recalling basic notions and results in connection with maximal 
regularity. In Section 3 we will review Mellin-Sobolev spaces, conically degenerate differential 
operators and their closed extensions. The above mentioned result on the existence of bounded 
imaginary powers for suitable extensions of the Laplacian on Mellin-Sobolev spaces will be estab¬ 
lished in Section 4. Section 5 contains an embedding result for the domains of the complex powers 
of the Laplacian. The subsequent section is devoted to the study of the porous medium equation 
in H®’^(B) and the proof of Theorem ll.il The principal result here is Theorem 16.II on the maximal 
regularity of the operators c — vA for c > 0 sufficiently large and v in the interpolation space 
Xi/q^q. Finally, in Section 7 we study the extension (11.61) in spaces with asymptotics and obtain 
Theorem 11.21 
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2. Preliminary Results on Parabolic Problems 

In this section let Xi ^ Xq be a densely injected Banach couple. 

Definition 2.1. For 9 G [ 0 , 7 r[ denote byV{9) the class of all closed densely defined linear operators 
A in Xq such that 

S'e = {z G C I I argzl < 0} U {0} C p(-A) and (1 + |z|)||(A + z)“^|| < Ag, z G Sg, 

for some Kg > 1. The elements in V{9) are called sectorial operators of angle 9. 

Given A G Vifi) we can define the complex powers A^ for Rez < 0 by a Dunford integral; 
composition with fc G N, then yields arbitrary complex powers, cf. Amann [51 III.4. 6 .5]. Of 
particular interest are the purely imaginary powers. 

Definition 2.2. Let A G V{9), 9 G [ 0 , 7 r[. We say that A has bounded imaginary powers if there 
exists some e > 0 and K > 1 such that 

G C{Xq) and ||A*‘|| < K for all t G [— e,e]- 

In this case, there exists a (j) > 0, the so-called power angle, such that ||A®*|| < for all t gM. 

with some M > 1, and we write A G 

The following statement is [22l Lemma 2.3]: 

Lemma 2.3. Let A G V{9) for some 9 > 0 and x G 'D{A‘^) for some 0 < (j) < 9. Then, for any 
Q <9' < 9 and 0 < p < fi, 

z i-A z'^A{A + z)~^x is bounded in Sg/. 

The proof actually shows slightly more, namely: 

Corollary 2.4. z >->■ z'^A{A + z)~^ is a uniformly bounded family in £(X>(A‘^), Aq). 

In UMD spaces (unconditionality of martingale differences property), the notion of the R- 
sectoriality, which is a boundedness condition stronger than the standard sectoriality, characterizes 
the property of maximal L^-regularity for the linear problem. 

Definition 2.5. Let 9 G [ 0 , 7 r[. An operator A G V{9) is called i?-sectorial of angle 9, if for any 
choice of Ai, ..., Xn G Sg, xi, ..., xn G Xq, N gN, we have 

N N 

II + ^^pIL2(o,1;Xo) - ^11 X!^P^I’IL2 (0,1 ;Xo)’ 

P=1 P=1 

for some constant C > 1, called the R-bound, and the sequence {cpj'^i of the Rademacher func¬ 
tions. Alternatively, we might have asked that, for some C >1, 

N N 

(2-2) II ^ + -^p) ^p|L2(04;Xo) — ^ II ^ ^I>^pIL2 (0.1 ;Xo)' 

P=1 P=1 

For later use, we state the following elementary observation. 


THE POROUS MEDIUM EQUATION ON MANIFOLDS WITH CONICAL SINGULARITIES 


5 


Lemma 2.6. Let A : Xi —>■ Xq be R-sectorial of angle 6, and let C be the R-bound. Denote 


S{9) = 


_ / sin(6») , e e [f ,7r) 


0 e[o,f) ■ 


Then, for any c > 0, A + c is again R-sectorial of angle 6 with R-bound < (7(1 + -c^)- 


Proof. For any Ai, Xn S Sg and xi ,xn G Xq, N gN, we have that 

N 

CpAp(7l + c + Ap) ^Xp 


p=i 


L2(0,1;A:o) 


N 


— ^ ^ £p(c + Ap)(A + c + Ap) 


p=i 

N 


L2(0,l;Xo) 


N 




CnXf 


p=l 


L2(0.1;Xo) 




1\ 


p^l 
N 

epXp 

p=i 


c + Ap 


N 

f I 'y ] £pc{A + c + Ap) 

p=i 

L2(0,l;Xo) 


L2(0,l;Xo) 


L2(0,l;Xo) 

where Kahane’s contraction principle has been applied in the last step, see [ITl Proposition 2.5]. □ 


Let A be a closed densely defined linear operator A : V{A) = Xi —>■ Xq. Assume that —A 
generates a bounded analytic semigroup on Xq. This is equivalent to the fact that c-\- A G V{6) 
for some c G C and some 6 > -k j^. Consider the abstract first order Cauchy problem 

. . / u'{t)+Au{t)= gif), tG{0,T) 

^ I u{0) = UQ 

in the Xo-valued L'^-space L‘^{0,T;Xq), where 1 < g < oo, T > 0. We say that A has maximal 
L'^-regularity, if for some q (and hence, by a result of Dore [TT], for all) the unique solution of (12.31) 
belongs to L^{G,T-,Xi) D VF^’‘?(0,T;Xq) D (7([0,T];Xi ) and depends continuously on the data 
g G L'^(0,T; Xq) and uq in the real interpolation space Xi := (Xi,Xo)i „. 

q ■Al q 

If the space Xq is UMD then the following result holds. 


Theorem 2.7. (Weis, [211 Theorem 4.2]) In a UMD Banach space any R-sectorial operator of 
angle greater than ^ has maximal L'^-regularity. 

For an alternative approach to maximal L'^-regularity based on a Hardy-Littlewood majorant 
type property of the resolvent of an operator instead of Rademacher boundedness, or for sectorial 
operators admitting a special type of operator-valued bounded i7°°-calculus, we refer to [22] and 
respectively. 

Remark 2.8. In a UMD space, an operator A G BIV{(j)) with (p < 1^12 is i?-sectorial with angle 
greater than 7r/2 by [71 Theorem 4] and hence has maximal L'^-regularity. This also is a classical 
result by Dore and Venni m- 


Next, we consider the abstract quasilinear parabolic problem of the form 

. . / u'{t) P A{u{t))u{t) = f{t,u{t)) + g{t), te(0,ro); 

^ I ir(0) = UQ 

in L^(Q,Tq-,Xq), such that V{A{u{t))) = Xi, 1 < q < oo, and Tq is finite. Maximal regularity for 
the solution of the linearized problem together with Lipschitz continuity will also imply existence 
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and maximal regularity for the solution of the original problem, as the following well known theorem 
guarantees. 

Theorem 2.9. (Clement and Li, [BJ Theorem 2.1]) Assume that there exists an open neighborhood 
U of uq in Xi such that A(mo) : Xq has maximal L‘>-regularity and that 

(HI) AeC^-iU,C{Xi,Xo)), 

(H2) /eCi-’i-([0,To] xC/,Xo), 

(H3) geL'}{0,To-,Xo). 

Then there exists a T > 0 and a unique u G L'?(0,T;Xi) n H()(0,T;Xo) fl C'([0,T];Xi ) solving 
the equation (12.41) on (0,T). 

3. Cone Differential Operators and Mellin-Sobolev Spaces 

3.1. Mellin-Sobolev spaces. For s G Nq, is the space of all functions u in Hpioc(^°) 

such that, near the boundary, 

(3.1) x^~~^{xd^ydy{u}{x,y)u{x,y))eLP(^[OA]xdM,^dyy j + |a| < s, 

where w is a cut-off function (near 9B), i.e. a smooth non-negative function w on B with a; = 1 
near 9B and w = 0 outside the collar neighborhood of the boundary. 

For an arbitrary s G M, dehne the map 

5., : ^ Cr(»"+'), v{t,y) ^ y). 

Let moreover kj : Uj C 9B —)■ R”, j = 1,..., iV, be a covering of 9B by coordinate charts and {q}j} 
a subordinate partition of unity. 

Definition 3.1. TLy^ (B)^ 5,7 G M, 1 <p < oo, is the space of all distributions on B° such that 

N 

(3.2) I|m||«;.^(b) = ^ ||5^(1 ® «;j)*(a;(/5jM)||/^.(Ri+n)-h 11(1-a;)u)||//|(B) 

i=i 

is defined and finite. 

Here, w is a (fixed) cut-off function and * refers to the push-forward of distributions. Up to 
equivalence of norms, this construction is independent of the choice of uj and the Kj. See [21 
Section2.2]. Clearly, all the spaces ?tp’^(B) are UMD spaces. 

We recall the following property for the Mellin-Sobolev spaces, which extends the standard 
Sobolev embedding theorem. 

Lemma 3.2. Let 1 < p < oo and s > {nl)/p. Then ||■u^;||^-,-y(B) < c||m||«j.7(b) ||u||.„=,(n+i)/ 2 (jj^ 
for suitable c > 0. In particular, ?tp"’'(B) is a Banach algebr^, whenever s > (n -|- l)/p and 
7 > (n-l-l)/2. Moreover, if s > {n-\-l)/p, a function u m?tp’^(B) is continuous onM°, and, near 
dE, 

|u(a:,y)|<cV-("+i)/2||u||„j,r(B) 

for a constant c' > 0. 

Proof. This is Corollary 2.8 and Corollary 2.9 in [23]. □ 

^Up to the choice of an equivalent norm. However, we will not distinguish the norms in the sequel. 
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The following is a slight improvement of Corollary 2.10 in [M]. We include the proof for the 
convenience of the reader. 

Corollary 3.3. Let 1 < p,q < oo, cr > 0, 7 S M. Then multiplication by an element m in 

n + l ^ + 1 

Hq ’ ’ ^ (B) defines a bounded map on for eaeh s G (—a,a). 

Proof. Let v G and denote by equivalence of norms. We can assume that v is 

supported near the boundary in a single coordinate neighborhood. Then 

||7n||^J + (n + l)/<!,(n + l)/2^j^ I ^ II Wj’'' (B) • 

Here, the first inequality is due to the fact that multiplication by functions in the Zygmund space 
defines a bounded operator in Hp for —tr < s < cr, cf. [501 Section 13, Proposition 9.10], and the 

second is a consequence of the fact that ^ Cf, cf. [501 Section 13, Proposition 

8.5]. □ 


Theorem 3.4. For s, 7 G M and 1 < p < 00 , the scalar product in yields an identification 

of the dual space to 'Hp’'^{M) with 'Hpfi’~'^(E), where p' is conjugate to p, i.e l/p+ 1/p' = 1. 

Proof. Follows by Definition 13.11 and the analogous result for the standard Sobolev spaces. □ 


The following result is [01 Lemma 5.4]. 

Lemma 3.5. Let sq, si, 70,71 G K, 0 < 0 < 1 and 1 < q < 00 . 






e,q ^ 


Wp'^-^{E) 

^s-5.7-e( 


Then, for arbitrary <5, e > 0, 
if q <2 

), ifq>‘^ 


with s = (1 — 0)so + Osi and 7 = (1 — 0)70 + O^i. 

Conversely we have, with the same notation and a very similar proof: 

Lemma 3.6. For arbitrary (5, e > 0, 

In the case of complex interpolation between Mellin-Sobolev spaces of the same weight, we have 
a sharp result on the order of the resulting Mellin-Sobolev space as follows. 


Lemma 3.7. Let sq) si, 7 G R and 0 < 6 < 1. Then, 

with s = (1 — 0)so + 0si. 

Proof. Recall that the statement is true for the standard Sobolev spaces in R" (see e.g. Equation 
4.(2.18) in [00]). Let a; = 1 on [0,1/2) and let {Vi,Ti)i with C B —)• R”, f = 1,..., M, be 
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a covering of B \ {[0,1/3) x 9B} by coordinate charts and {ipi} a subordinate partition of unity. 
With the notation as in Definition o we have the following equivalence 


u G ■H;’^(B) 

<5=^ Sj{l ^ (ri)*(^i(l - a;)u) G Vi,/ 

^ S^il^KjUoj^ju), (t,)*(V'*(1-w)w) G Vi,j 

^ ecu, (l-u;)uG [i«;°’^(B),H;^'^(B)]e 
^ uG [//«“’^(B),7/;-'^(B)],, 

where we have used the linearity of the push-forward together with Proposition 4.2.1 in m- □ 


3.2. Cone differential operators. A cone differential operator or conically degenerate operator 
A of order ^ is a differential operator of order /r on the interior B° of B which, in local coordinates 
near the boundary, can be written in the form 

(3.3) A = x~^'^^aj{x){—xdxy with a^-G (^“([0,1), Diff^”-^(5B)). 

j=o 

The operator A is called B-elliptic (or degenerate elliptic), if the principal pseudodifferential symbol 
cr^iA) is invertible on T*B° and, moreover, in local coordinates {x,y) near the boundary and 
corresponding covariables (^, ??), the rescaled principal symbol 


(3.4) 


x^(^'^iA)ix,y,^/x,y) 




1=0 


is invertible up to a; = 0. 

Moreover, one associates to a cone differential operator its conormal symbol. This is the operator 
family ctm(A) : C ^ defined by 

(3.5) aM{A){z) = aj (0)zT 

1=0 


One is mostly interested in the points where aM{A) is not invertible. For this, the precise choice of 
s and p in (13.511 is irrelevant due to the spectral invariance of differential operators in these spaces; 
it is often convenient to work with s = 0 and p = 2. 

The Laplacian A induced by the metric g is a second order cone differential operator. Near the 
boundary A can be written in the form 

(3.6) A=-^{{xd^y+ {n-l){xdx) +^d), 

where Ag is the boundary Laplacian induced by h. Hence A is B-elliptic. The conormal symbol 
(Jm{A) of A is given by 

(3.7) (Jm[A){z) = - [n-l)z + Ag. 


A conically degenerate operator A acts in a natural way on the scales of weighted Mellin-Sobolev 
spaces: 

is bounded for all s, 7 G M, 1 < p < oo. 
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3.3. Closed extensions. Sometimes it is necessary to consider a cone differential operator A as 
an unbounded operator in a fixed space If A is IB-elliptic, then the domain of the minimal 

extension, i.e. the closure of A considered as an operator on C^(B°), is 

iJ- 

(3.8) e n G 

e>0 j—0 

If, in addition, the conormal symbol of A is invertible for all z with Re z = (n -I- l)/2 — 7 — /r, then 

The domain of the maximal extension, defined by 

2?(^s.max) = {UG : Au G 1W*’^(B)}, 

satisfies 

T’(As^niax) = T>(As_i„in) © 

where f is a finite-dimensional space consisting of linear combinations of functions of the form 
x~P log^ X c{y) with p G C, fc G No and a smooth function c on the cross-section. The space £ can 
be chosen independent of s. This result has a long history, see e.g. 0, US], ESI, EZI; the present 
version is due to Gil, Krainer and Mendoza [13] . 


3.4. Extensions of the Laplacian. We are interested in the values of z where (Tm(A) is not 
invertible. We denote by 0 = Aq > Ai > ... the distinct eigenvalues of Ag and by Eq, Ei,... the 
corresponding eigenspaces. Moreover, we let iTj G £(L^(i9B)) be the orthogonal projection onto 
Ej- 

The non-bijectivity points of aM{^) are the points z = and z = with 



Note the symmetry q'^ = {n — 1) — q^ . It is straightforward to see that 


(3.10) 


(z^ - (n - I)z -I- Ag) 


-1 


00 

E 




j=o 


Hence, in case dimB 7^ 2, where the q^ are all different, the inverse to (Tm(A) has only simple 
poles in the points q^. For dimB = 2 the poles at qf,j^ 0, are simple, while there is a double 
pole at qQ = q^ = 0. 

With q^, j ^ 0, we associate the function spaces 

£ ± = ujx~'^^ ® Ej = {w(a;) x~‘^i e{y) : e G Ej}, j G N. 

For j = 0 we let 


(3.11) 




(jj ® Eq + a; logoi © Eq, dimB = 2 
ujx'^o © Eq, dimB 7^ 2 


For later use note that A maps the spaces £ ± to (^((“(B”). 

^3 

Furthermore, we introduce the sets /.y, 7 G K, by 

/^ = j G No} n - 7 - 2, ^ - 7[. 

The following is Proposition 5.1 in m combined with Theorem 3.6 in m- 
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Proposition 3.8. The domain of the maximal extension of A in is 

= 2^(As,min) © ^9*’ 

qfely 

In case — 7 — 2 for all j, the minimal domain is 2?(As_i„in) = {^) ■ 

Corollary 3.9. The domains of the closed extensions of A are the sets of the form X>(As^min) ® f, 
where £ is any subspace of 0 ± _, £ ±. 

9j <^-‘1 Hj 

Definition 3.10. Given a subspace £± of £ ±, we define the spaee f'i as follows: 

i) If either qf ^ 0 or dim(B) ^ 2, there exists a unique subspace C Ej such that £± = 

^ J J Hj 

® E_y Then we set 

£\ =ujx~'^'^ 

'ij •’ 

where is the orthogonal complement of in Ej with respect to the L^{dM)-scalar 
product. 

ii) For dim(B) = 2 and q^ = 0 define = {0} if ^ = £ 0 , ^o" = £0 if ^0 = {0}, and 
^o" = ^0 */fo = ^ © ^0- 

Note that £\ is a subspace of £ ^. 

We will now confine ourselves to closed extensions of A in (B) with domains 
P(AJ = Il(A,.^i„) © © f„± c 

g±G7^ ' 

chosen according to the following rules: 

(i) If qf eljC] I-^, then 

(ii) If 7 > 0 and q^ € Ij\ then £^± = £^±. 

(hi) If 7 < 0 and qf G I.y\ then £± = {0}. 

J Hj 

In particular, 21(A J = 2?(As,i„ax) if 7 > 1 and I1(AJ = X>(As,min) if 7 < -1. 

Theorem 3.11. Let s >0, 6 G [ 0 , 7 r[, and (j)> 0. For I 7 I < dim(B)/2 let A^ be an extension with 
domain in chosen as above. Then c — A^ G 'P{0) (7 BTV{4>) for suitably large c > 0. 

Proof. For s = 0 this follows from Theorem 2.9 and Remark 2.10] with Theorems 5.7 and 4.3 
in [5^. For s > 0 we apply Theorem 3.3 in [M]. □ 

4. Bounded Imaginary Powers for the Laplacian on Negative Order 

Mellin-Sobolev Spaces 

Let X be a complex Banach space and X* its dual. The adjoint A* of a sectorial operator 
A with bounded imaginary powers will also be sectorial on the same sector and have bounded 
imaginary powers with the same power angle, provided T>{A*) is dense in X*, see Propositions 
1.3(v) and 2 . 6 (v) in [TU]. We will use this fact in order to extend Theorem 13.Ill above to the case 
of Mellin-Sobolev spaces of negative order. 
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Theorem 4.1. For s > 0, I7I < and 1 < p < 00 consider the closed extension A_g of the 


Laplacian with domain 
(4.1) 


25(A_ J = 0 0 




where the asymptotics spaces £± satisfy the conditions (i), (ii) and (iii) in Definition I3.10L Then, 
for any 9 € [0, tt) and (j) > 0, there exists some c > 0 such that c — A_g € V{9) n BI'P{(j)). 

Proof. We will show that (ED is the adjoint of a closed extension of the Laplacian on another 
Mellin-Sobolev space, namely, of the extension p, with the domain 


(4.2) 


~ T^{As,-~i ,p' ,m\n) ffi —0*’ 




where l/p+ 1/p' = 1. Here X>(As^_..),_p/^min) denotes the minimal domain of A in T-Lp, "’^(B). 

In fact, by Theorem 5.3 in m. the domain of the adjoint of is 

Il((A,,_^,p,)*)=^(A-..™i„)© 0 

qfel-, 

This coincides with (14.11) . since, by the definition of the orthogonal complements, 

(4)^=£,*- 

In view of the above stated result on the adjoints of sectorial operators with bounded imaginary 
powers, it suffices to show that, given any 9 G [0, tt) and any ^ > 0, the operator c — Ag^-^y with 
domain (|4.2I) is sectorial of angle 9 and has bounded imaginary powers with angle cj), provided c 
is large enough. This in turn will follow directly from Theorem 13.111 provided the asymptotics 
spaces satisfy the conditions (i), (ii), and (iii) in Definition 13.101 

In order to see this we note first that + q~ = n — 1, so that q^ G I.y if and only if qj G I-^. 
As a consequence, we may rewrite the direct sum on the right hand side of (14.2|) as the direct sum 
over all qJ G I-^. It remains to check that the asymptotics spaces have the properties in (i), (ii) 
and (iii). We distinguish the three cases 

• We have qf G D /_.y if and only if qf £ Cl I-j. Hence the assumption £\ = £j for 

J J Qj Qj 

A_s shows that also = £± , so that (i) holds for the extension with domain ( 03 ). 

Qj Qj 

• If 7 > 0 and qf G then qf G I.y\I-.y. Since A_j, satisfies condition (ii), £ ± = £ +. 

J J Qj Qj 

Therefore £"*7 = {0}, which says that the domain (14.21) satisfies (iii). 

Qj 

• Similarly, if 7 < 0 and qf G /_-y \ , then qf G I.y\ . As (iii) holds for A _^, we have 

£ ± = {0}. Therefore = £ ±, and the domain (14.2p satisfies (ii). 

Qj Qj Qj 

□ 


4.1. Application. By assumption, Ai is the largest negative eigenvalue of the Laplacian Ag in¬ 
duced by the metric h on the boundary. We let 


_ n — 1 
=--1 


n — 1 


)2 - Ai = > 0. 


(4.3) 


2 
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We choose the weight 7 in the interval 

tAA\ n - 3 

(4.4) —-— < 7 < mm 


n — 3 
2 


_ n + 11 
En, 2 / 


and consider the closed extension of the Laplacian in the Mellin-Sobolev space := Xq, 

1 < p < 00 , s G K., given by 

(4.5) P(A J = H;+2’^+2(1) © C := Ai. 

Note that ^^— 7—2 = ^^—7 G (—£n, 0) does not coincide with one of the qf so that 

is the minimal domain of A in 'H®’'>'(B). In addition, the choice of the single asymptotics space C, 

the constant functions, meets the requirement in Theorem 13.Ill 

We can therefore deduce from Theorem 14.11 that for any cj) > 0, the operator c—A^ with domain 
(B31) and 7 satisfying (14.4p has BXV{(j)) provided c is sufficiently large. But more is true. Let 
c ^ IR<o- Based on the results in it was shown in [531 Theorem 4.1] that c—A^ : H®+^’'’''''^(B)©C 
is invertible for s > 0. Since it follows from Corollaries 3.3 and 3.5 in [5S] that the kernel and the 
cokernel are independent of s and p, we find that it is also invertible for s < 0 . 


Theorem 4.2. For any c > 0, 0 G [0, tt) and ^ > 0, the operator c — with domain (14.51) and 7 
satisfying (14.4|) belongs to V{9) n 


Proof. As pointed out above, we know the result already for large c. Since the resolvent exists 
outside K<o and 9 G [ 0 , 7 r) is arbitrary, we have c — A^ G 'P{9). According to Lemma IIL4.7.4 in 
[2] we find for any </) > 0 a c' > 0 and M > 0 such that the norm of (c' — A^)^ in is 

< for —1 < Rez < 0. Denote by Lg the positively oriented boundary of the sector Sg. 

Using the Dunford integral for the complex powers, see III.4. 6 .5 in [5], we see that 


ll(c-AJ^||<||(c' 
< Ale'll ^“^1 + 


-AJ^II + ll(c-AJ^ 


27r 


[ (-A)^(c- 


-{c-A^rw 

A, + A)-1(c' 




A, + X)-^dX 


for a suitable constant M'. Since 9 can be taken arbitrarily close to tt, we have c — A^ G BXV{(j)) 
by Lemma III.4.7.4 in [2]. □ 


5. Embeddings for the Domain of the Complex Powers of the Laplacian 


At this point, we are able show how the interpolation between a Mellin-Sobolev space on one 
hand and a direct sum of a Mellin-Sobolev space and the constant functions on the other can be 
controlled. As a consequence, standard theory for sectorial operators furnishes some information 
about the domain of the complex powers of the Laplacian. We first prove the following. 

Lemma 5.1. Let 5,7 G M and p G (l,oo). Suppose u G (B) n and that in local 

coordinates near the boundary, xdxU and dy.u belong to 77®+^’'’'''"^(B). Then u G 77(As_niin) + C. 

Proof. From (O and the assumption we conclude that Au G 77®’'’'(B). Therefore, u belongs to 
the maximal domain of A in 77®’'’'(B), and Proposition 13.81 shows that 

U G 77(As^inin) ® ^9*' 

From the above direct sum and the assumption on dxU we see that only the asymptotic component 
that corresponds to = 0 can occur, which completes the proof. □ 












THE POROUS MEDIUM EQUATION ON MANIFOLDS WITH CONICAL SINGULARITIES 


13 


We can use the above lemma to control certain interpolation spaces as follows. 

Lemma 5.2. Let s G K, p, g S (1, oo) and 6 S (0,1). Then, the following embeddings hold 

+ c -A + c)g,g ^ + c, 

for every e > 0. 

Proof. Concerning the first embedding we know from Lemma 13.61 that 

By standard properties of interpolation spaces (see e.g. Proposition 4.2.1 in [H]), we also have 
(H*’^(B),7^;+2''''+2(B))e., + c)g,g. 


Therefore, 

^s+2e+E.7+2e+e(B) ^ (?^;-t'(B),h;+2>t'+ 2(B) +C)e,,. 

Moreover, 

C -A (H;’'>'(B),7^;+2'T'+2(B) +C)e.,. 

Hence, we conclude that 


^s+2fl+e.7+20+e(B) + C -A (^^’^(B), ■H;+2’'>'+2(B) + C)e,,. 

Let us show now the second embedding. For any 

u e (h;’'>'(b),h;+2'T'+2(b) + c)e.„ 

xdxU and dy^^u belong to 

for any e > 0 by Lemma 1531 Then Lemma lOI and (13. 8|) show the assertion. □ 

In this way, we gain some information about the domains of the complex powers of given in 

(14.5F From [2j 1.(2.9.6)], [2j 1.(2.5.2)] and the previous lemma, we obtain the following. 

Corollary 5.3. Let s G M., j as in (14.4F p G (l,oo), 0 < Rez < 1 and c > 0. For the closed 
extension A^, given in (HAl) . we then have 

^,+2Rez+e.7+2Re^+e(B) ^ ^ ^,+2 Re z-e.7+2 Re z-e (g) ^ ^ > q_ 

6. The Porous Medium Equation on Manifolds with Cones 

We will treat the porous medium equation as a quasilinear problem. Using the multiplier 
properties of the Mellin-Sobolev spaces of sufficiently large order and weight, we will apply the 
maximal regularity theorem by Clement and Li. We denote by (g*-^) = {gij)~^ and (h*-^) = 
the inverses of the metric tensors of g and h in local coordinates. Then the following identity holds 

A(u"*) = mu'^~^Au + m(m — 1 )u™“^(Vm, Vu)g, 


where 


Vu = > g 


du d 


^ dx'^ dxt ' 


and, in local coordinates (x, y) near the boundary, 

{Vu,Vv)g = ■^(^{xdxu){xdxv) 


E 


dy^ dyi J 
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Equations (fTTl) . (ra then take the quasilinear form 

(6.1) u'{t) + mu"^~^{t)Au{t) = /(u, <) — m(m — l)u"*“^(t)(Vu(t), Vu(t))g 

(6.2) u(0) = uq- 

In the following, we will show maximal LP-regularity for the linearized term of the problem 
for appropriate initial data. The proof is inspired by Theorem 5.7 in m- We first construct a 
parametrix with the help of a suitable partition of the space and a Neumann series argument. 
Then we apply i?-sectoriality perturbation arguments. 


Theorem 6.1. Let A^ be chosen as in (1431) and the weight 7 as in g31). We assume that s > 0 and 
that p and q are so large that _|_ 2 ^ and 7 > + - . Then for any u S Xi := (Xi, Xq) 1 

such that u > a > 0 onM and sufficiently large c > 0 , c — uA^ : Xi —>■ Xq is a well defined closed 
linear operator that is R-sectorial of angle 6 for any 9 in [0,7r). 

The statement extends to the case where — 1 + + | < s < 0 with p, q as before. 

Proof. Recall that Xq = Tlp’'^{M) and Xi = 0 C. By Lemma lOl 

(6.3) Xi „ := (Xi,Xo)i „ -o ’""(B) © C 

q'^ q'^ 

for any e: > 0. In the sequel, we will only use this property of u. 

Since s + 2 — | > |s| + and 7 + 2 — | > 21^, Lemma [!0 

implies that u is continuous and constant on i9B; according to Corollarv l3.31 multiplication by u 
induces a bounded map in 'H®’'^(B). Thus, uA^ is a well defined closed linear operator in Xq with 
domain Xi. 

For z G B, let u{z)A^ be the realization of u{z)A in Xq with domain Xi. For every c' > 0, 
c' — Aj is i?-sectorial of angle 9, for each 9 G [0, tt) by Theorem l4.2l and [TJ Theorem 5]. Fix c, c' > 0 
with c' < c/||m||oo- Denote by Cp the p-th Rademacher function. The i?-boundedness of c' — A^, 
cf. (12.111 . (12.211 . implies that for any Xi,...,Xk G {2 G C | |argz| < 9}, K G N, and xi,...,xk G Xq, 


K 


Xp{c - u{z)A^ + Xp) ^epXf, 


p=\ 


L2(Q,1;Xo) 


+ ^-^ + c') (c'-A, + ^-c' + ^) 
u{z\ u(z) /V u(z) u(z)J 


l{z) 


— c 


< 


n 

p-i 
K 

^u(z) ' u(z) 

P—1 ^ ^ ' 


Xp 


u{z) u{z) 
X 




u[z) 

Xn \-l 


l{z)- 


ioX, 


L2(0,l;Ao) 


u[z) u{z) 


) 


^p'^ P 


L2(0,l;Xo) 




.u{z) 


P=1 


< 


A 

{c\c') + (^ - c')||(c' - AJ-i|U(^„)C"(c')) I 

'' ' p—1 


u{z) u{z)t 

K 

CnXn 

L2(0,l;Xo) 


eoXc 


L2(0,l;Xo) 


for constants C"(c'), C'"(c') > 0 that depend only on c'. Hence, for any c > 0, c — u{z)A^ is 
i?-sectorial of angle 9 and the i?-bound is uniform with respect to z. By a well-known perturbation 
result, see Theorem I in [TB], the sum c — u{z)A^ — v(z)A^ is again i?-sectorial of angle 9 with a 
slightly larger i?-bound uniformly in z, provided the norm of v in Xi is sufficiently small, say 
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ll'f^llxi < Here we use the fact that multiplication by v induces a bounded map in Xq, by (16.3|) 

q 

and Lemma o 

In order to establish the i?-sectoriality of c — uA^ assume first that s = 0. For r > 0 choose 
an open cover of B, consisting of balls Bj = Br{zj), Zj G B°, j G {1, N}, not intersecting the 
boundary, together with a collar neighborhood Bq = [0,r) x 9B. We may assume that also the 
balls B 2 rl 2 {zj) do not intersect 9B. Let w : R [0,1] be a smooth non-increasing function that 
equals 1 on [0,1/2] and 0 on [3/4, oo). Denote by d = d{z,z') the geodesic distance between two 
points 2, z' G B with respect to the metric g. We define 

Uj{z) = w u(z)-h - w u(2j), J = and 

d{z,dM)\\ 

and write 


uo{z) = u;(^^^^)u{z)+(l-w( 


c - Uj{z)A^ = c - u{zj)A^ + {u{zj) - Uj{z))A„ 


where each realization is considered with the same domain Xi as A^. Since ||u(2j) — itj(2)||c(B)) 
and therefore the norm of u{zj) — Uj{z) as a multiplier on can be made arbitrarily small by 

taking r small (for j = 0 recall that u is constant along the boundary), each c — UjA^ is i?-sectorial 
of angle 0. 

Fix a partition of unity (j)j G (^““(B), j = 0, subordinate to the Bj and functions 

ipj G C°° (B), supported in Bj , with = 1 on supp (j)j. 

Let / G Ai, 5 G Ao and A G C. Multiplying 


A/ - uAJ = g 


by (j)j "we find 


Mjf - u^{(j)jf) = (j)jg - [uA,, 


Applying the resolvent of Uj A^ we get 

= (A - UjAJ-Vi5 - (A - UjAJ-i[uA^, 4>j]f. 


In view of the choice of ipj, the above equation becomes 

hf = V’j(A - UjAJ"Vj5 - '</’i(A - MjAJ"^[uA^, 

Summing up we obtain 

N N 

(6.5) / = X! - X! “ UjAJ-^[uA^, f/j]/. 

j=o j=o 

The commutator [uAs.,4’j] is a first order cone differential differential operator with non-smooth 
coefficients. In view of (16.311 and the subsequent considerations, it maps D(A^) to {V>) for 

suitably small d > 0 (note that (j>o is constant near the boundary). According to Corollarv l5.31 the 
latter space embeds into T>{{c — A^Y) for any v < 5/2. We may now apply Corollary 12.41 and find 
that 


(6.6) (A - Uj AJ ^[uAg, (/j] 0, A ^ oo, 

in £(Ai). Hence, A — uAg will be left invertible whenever A ^ IR<o is sufficiently large. 







16 


NIKOLAOS ROIDOS AND ELMAR SCHROHE 


We denote the left inverse by Sx. Then 

N 

(A - = (A - uAJ A - [mA^, (I)j]Sx) 

j=o 

N N 

= - [uA^^'/'jIS'a) - y][uA3,V’j](A - - [mA^,(/)j]S'a). 

j=0 j=0 

Since ~ ^ E^qI'^As; </>j] = 0, we conclude that, on Xg, 

N 

(6.7) (A-wAJS'a = I - UjAs)~^{<Pj - [wA,,<))j]S'a). 

j=o 

We will next argue that the sum on the right hand side tends to zero in C{Xg) as A —)> cx). Choose 
1/2 < z/' < 1 and write 


(A-UjAJ i = (c-MjAJ '^'(c-UjAJ(A-MjAJ ^(c-m^AJ . 

Then the product of the last three terms tends to zero in C[Xg) as A —>■ oo by Corollary 12.41 The 
operator (c — UiA^)~''' maps Xg into T>((c — UiA^)''') ^ 0 C by Corollary 15.31 As 

the commutator [uA^,ijjj] is a cone differential operator of order one it maps 0 C 

continuously to Xg. Hence the sum on the right hand side of (16.71) tends to zero and thus A — uA^ 
has also a right inverse for large A ^ K<o. Let 

N N 

-R(A) = and Q(A) = Yi’ji>^ “ UjA^)~^[uA^,cl)j]. 

j=o j=o 

As Q(A) —0 by (16.61) . we conclude from (16.51) via a Neumann series argument that, for A large, 

OO 

(6.8) (A - uAJ-i = ^(-l)'=Q'=(A)i?(A). 

/t-0 

Starting from dSSl) and splitting off the term for k = 0 we have 


K 


(6.9) I y] Ap(Ap 0 c-mAJ Lpe^ 

p=i 

K N 

- I Xpei_ 

p=l j=g 

OO K N 

^ ^ I ^ y ApQ (Ap + c) ^ '0j(Ap + c — Uj Ag) (j)^Xpe^ 


L2(0.1;Xo) 


L2(0,l;Xo) 


k=l p=l 


j=0 


L2(o,l;Xo) 
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We then estimate the first term on the right hand side by 

N K 

- X! X! “ UjAJ-Vja;pep 


i=o p=i 
N 


L2(0,l;Xo) 


K 


Xn€r 


j=0 


p=l 


L2(0,l;Xo) 


K 


< Ci{N + 1) max || 'Y, 4‘jXpe^ 


p=i 


l2(o,1;a:o) 


K 


( 6 . 10 ) 


< C2(iV + 1) I ( XpCf, 

P=1 


L^{0,l-,Xo) 


as required. Next we focus on the summands in the second term for fc = 1, 2,.... They are 

N K 

(6.11) < I ApQ (Ap + c)'0j(Ap + c — itj Ag) (f>jXp€p 

j=0 " p=l LH0,1-,Xo) 

Expanding the product, we regroup this into a sum of k factors; the first equals 

(6.12) ^jAp(Ap + c - MjA^)“^ 
the others are of the form 

[uA^, (j)j]ijji{Xp + c - UiA^)~^ = [uA^,(j)j]4!i{c - UiAJ“^(c - UiA^){Xp + c - m^A^)”^ 

for suitable i,j G {0,...,A}. Now [uA^, cj)j]tpi is a first order cone differential operator. By a 
similar argument as before, the norm of [uA^, (pj]ipi(c — UiA^)~^ in C{Xq) will tend to zero as 
c —oo. Given e > 0, we choose c so large that all these norms are < e. Using successively the 
i?-sectoriality of c — ujA^ together with Lemma 12.61 and (12.211 we can estimate (16.111) by 

K 

< CsiN + 1){{N + l)Ce)'^-i ^ XpEp 


p=i 


L2(0,l;Xo) 


with a hxed constant C. Summing over k and combining this with the estimate (Ib.lOL we see that 
(lOll is 

K 


<C5 


A+ 1 


1 - (A + 1)C£ 


1^^ 

p=i 


•p^p 


L2(0,l;Xo) 


which establishes the i?-sectoriality in Se- 

Next we will show that the resolvent of c — uA^ for s > 0 is the restriction of the resolvent of 
c — itAg to It suffices to prove that the following identities 

(uAg — A)“^(uA — X) = I and (uA — A)(uAg — A)“^ = I 

hold in Xi and Xq respectively. Let x € (D(uAq) = 0 C such that uAx G 'H®’'’'(B). 

I 0C, for all £ > 0, which is a Banach algebra. 


s~l“2___ s '*Y“I“2___ £ 

By (16.(111 . Xi embeds into T-Lp ’ ’ " ( 

a 


,s+2-^-e.7+2-^-£ 


As shown in Lemma below, also u ^ belongs to Hp 
elements in that space are multipliers in 'Hp'^{ 

X G n (B) 0 C), and hence x G 'Hp+^’'’'+^(B) 0 C. 


we conclude that Ax G H®’'’' 


, for all £ > 0. Since 
Therefore, 
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Let US now treat the case s > 1, s S N, by induction. Let uq, ...,yn be local coordinates on 
the support of (j)i, with the understanding that, in the collar neighborhood of the boundary, we 
use local coordinates x,yi, ...,yn and replace the derivative dyg in the computation below by xdx- 
Denote by M^. the multiplication operator by 4>i. Then 


K 


'y ^pi^p “b ^As+i) ^p^f 


p=i 


K 




^Ap(Ap + c-uAJ ^xpc^ 
p=i 


L2(0,l;^p 


N K 


— ^ ' I y ' dyM(f,.Xp(\p + c — uAg) Xp€p 

|a|<l i=0 p=l 

N K 

|a|<l j=0 p=l 

N K 




L2(0,(B)) 


|a|<l 3=0 p=l 


L2(0,1;WJ'''(B)) 


K N 


L2(0,(B)) 


~ 'y ^ I ^ ^pi^p ~b dyM^.XpCp 

|a|<l P=lj=0 
K N 

+ Y. IX] X] ^p^^p + ^ “ 

|a|<l p=li=0 

(6.13) X [d‘^M^.,c - uAJ(c - uAJ-^(c - uAJ(Ap + c - uAJ-ia;pe^ 


L2(0.1;WJ''^(B)) 


By the assumption on q and p, the coefficients of the second order cone differential operator 
[dyM^^.c — uAg] are multipliers in '^*’'''(8) (note that by the induction hypothesis now u belongs 

s+3—— —£ 

to Up ’ ’ ^ (B) © C). Therefore, the operators [9“M0.,c — uA^]{c — uAj,)“^ are bounded 

maps in ^*’'’'(18). Using the i?-sectoriality of c — uA^ in Hp’^(B), we obtain from (16.131) 

K 

y ' •^p(^P + C — MAg_|_]^) XpCp 


p=i 


L2(0,1;-HJ+LT(b)) 


K N 


K N 


— ^6 X I X X dyM^.Xptp 

|a|<l p=li=0 




■^6 X ^ 

|o|<l P=1 j=0 


x(c - uA^)(Ap + c - uA^) ^XpCp 


L2(0,1;WJ’''(B)) 


K 




•p^p 


P=1 


+ Ce max 

L2(o, 1;'HJ+1'T'(B)) i,a 


[d^M ^., c - uAJ(c - uAJ- 


K 




P=1 


Xot, 


L2(0.1;WJ’'^(B)) 
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< 


< 


K 

Cq I ^ ^ 

p=l 

K 

II ^ ^ ^p^p 

P=1 




K 

+ CtI y^a;pep 

p=i 






for suitable constants Ce, Cr and Cg. 

The case where s € K, s > 0, follows by interpolation, see Lemma iTTl and [151 Theorem 3.19]. 

Assume finally that —1 + + | < s < 0. We basically proceed as in the case s = 0. 

The crucial step is to prove that the norm of u{zj) — Uj(z) as a multiplier in can be 

made arbitrarily small by choosing the radius r of the covering small. In view of the fact that 
|s| < 1 — (n+l)/p, Corollary 13.31 implies that it is sufficient to show that the norm of u{zj) — Uj{z) 
in T-Lp tends to zero as r —>■ 0 for any fixed £ > 0 (note that the constant part of 

u cancels when we take the difference). In view of the fact that our assumptions on s,p and q 
imply that u € for suitably small (5 > 0, this will follow form a standard 

interpolation inequality, see e.g. Proposition 1.2.2.1 in [5], provided we show that 


(6.14) \\u{Zj) - Mj(z)||^0,(n + 1)/2(JJJ 

(6.15) \\u{zj) - Mj(z)||^i,(„+i)/2(jj) 

Now (16.141) is obvious. As for (16.151) write 

(6.16) - 


tends to zero as r 0; 
is bounded as r ^ 0. 


Uj{z) - u{zj) = “ u{Zj)). 

Being an element of for some d > 0, m is Lipschitz continuous in the 

interior of B, so that 

\u{z)-u{zj)\ < Ld{z,zj) 

for some L > 0. In view of the fact that ^ is bounded uniformly in r, we obtain a uniform 

bound on the derivatives in (16.161) on the balls in the interior. Concerning Bq, we notice that, for 
z = {x,y) we can take d(z,9B) = x in the definition of uq and use the fact that dyUJ^^J = 0, 

while ^ is again bounded. This completes the argument. □ 


Lemma 6.2. Let u € Hp‘’’^“(B) 0 C, for some sq > 7 o > If u pointwise invertible, 

then u~^ G 'Hp°’'^°{M) 0 C. As a consequence, 'Hp”’’^”(B) 0 C is spectrally invariant in C'(B) and 
therefore closed under holomorphic functional calculus. 

Proof. The case, where sq is a positive integer, is straightforward. Otherwise we choose the Bj, uj 
and (pj as in the proof of Theorem 16.11 We assume that r is chosen so small that 

\\Uj{z) - u{Zj)\\sup < ^\u{Zj)\. 

In an initial step, we will show that the inverses of the Uj have the asserted property. We 
consider first an interior ball Bj, j > 1. We may assume that even B^j./ 2 {zj) is contained in single 
coordinate neighborhood for B°. The function Vj{z} = Uj[z) — u{zj) is supported there, and its 
push-forward under the coordinate chart belongs to iLp“(]R"+^). To simplify the computation, 
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below, suppose that u{zj) = 1. Theorems 6 and 10 in [3] show that Vj{l + Vj) ^ also belongs to 
Its pullback under the coordinate map then is an element of Since 

Uj^ = (1 + Vj)-^ = 1 - + Vj)-^ 

we conclude that uj^ belongs to © C. 

As for Vo, we first note that it also is an element of Using possibly a further 

partition of unity, we may assume that it is supported in a single boundary chart. Denote by u* 
its push-forward under the coordinate map. According to Definition 13.11 V = S(n+i)/ 2 '<^* is an 
element of By the above theorems, the same is true for U(1 + V)~^, and we deduce 

that z;o(l + uo)“^ is an element of Now we write 

■ao ^ = (1 + = 1 - uo -f uo(uo(l -f uo)“^). 

Lemma 132] then shows that the right hand side is in © C. 

To see that u~^ also has the asserted property, we recall that Uj coincides with u on the support 
of (pj, so that 

(6.17) 1 = (pjuj^u. 

This shows the spectral invariance of ?^p°’'’'°(B) © C in C'(B). 

It is well-known that this implies the closedness under holomorphic functional calculus. We recall 
the argument for the convenience of the reader. If h is a holomorphic function in a neighborhood 
of Ran(—u), then 

^(^) = 7:“ f h(—X)(u + X)~^d\, 

2m Jr 

where T is a finite path around Ran(—u). Since the spectral invariance implies the continuity of 
inversion, we conclude that h{u) G 77®”’'’'“(B) © C. □ 

Lemma 6.3. Under the conditions of Lemma \6.‘A let U be a bounded open subset o/?7®“’'’'“(B) ©C 
consisting of functions u such that Reu > a > 0 for some fixed a. Then the subset {u~^ \ u G U} 
o/'H®“’’'’'“(B) © C JS also bounded. 

Proof. If So is a positive integer, this is straightforward from the identity dy{u~^) = u~^dyU. For 
non-integer sq we transfer the problem to as in the proof of Lemma l6.2l and apply [U Theorem 
11 ]. □ 


Remark 6.4. Let G K and 6 in [0, tt). Under the conditions of Tfeeorem l6.ll the operator c—u'^A^ 
is a closed linear operator on Xq with domain Xi which is R-sectorial of angle 9, provided c > 0 
is sufficiently large. 


g _j_ 2 _—_^ ”|”2___ s 

In fact, since u belongs to Up ® (B) © C for each e > 0, the same is true for u'^ 

by Temma l6.2l As pointed out in the beginning of the proof of Theorem \6. 11 this is all we need to 
know about u for the proof of the R-sectoriality. 


Summarizing what we have found we obtain the following result for the porous medium equation. 


Theorem 6.5. Let Ai be the largest nonzero eigenvalue of the boundary Laplacian Aq, induced by 
the metric h on i9B. Recall that dim(B) = n + 1, that 


^71 


n — 1 



> 0 , 


2 
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and that 7 satisfies (03), i.e. 


n — 3 

—-— < 7 < min 


n — 3 
2 



Next choose p, q so large that 

(6.18) + and 7>^ 

p q 2 


4 


Then, for any s > — 1 + | and for any strictly positive initial value 

there exists some 0 < T < Tq such that the porous medium equation dll), dH) in the space 
L^(0, T;'H®’'>'(B)) has a unique solution 


Proof. We shall apply Clement and Li’s Theorem l2.91 to the porous medium equation in the form 
(16.11) . (16.21) . with Xq = and Xi = Hp+^’'’'+^(B) © C. 

For each initial value uq € Xi/g g = (Xi, Xo)i/g^g such that mq > a > 0, the operator 
has maximal L^-regularity by Bemark |6.4l It remains to check the conditions (HI) and (H2). By 
our assumptions on s,p, q and 7, (lOl) in connection with Corollary 13.31 implies that Xi/g g embeds 
into a Banach algebra of bounded multipliers on Xq. So let [/ be a bounded open neighborhood 
of Mo in Since {Xi, XQ)ifg g embeds into C'(B), we may assume U to consist of functions 

with real part > a. For any mi, U 2 G 1/ and real v, we have then 

(6.19) M^ - M 2 = (m 2 - Mi);^ /(-A)^(mi + A)“^(m 2 + A)“^iiA, 

2711 Jr 

where F is a fixed finite path around U„g£/Hnn(— m) in {Re A < 0}. 

Concerning (HI): Equation (16.191) in connection with the embedding (16.31) . Lemma Lemma 

El 

and the boundedness of U implies that 

\\muT-^A-mu^-^A\\cix,,Xo) < \\muT-^ - mu^^WciXo) 

(6.20) < Ci\\ui — U 2 \\c{Xo) ^ (^ 2 \\ui — U 2 \\xi , 

q 

for suitable constants Ci, C2 > 0. 

Concerning (H2): Equation (16.31) shows that Vmi and Vm 2 belong to TLp^^ 2/g 

for each £>0. As2/(jf<l, they are elements of Xq. Write 

||/(Mi,ti) -m(m- 1 )m™“^{Vmi, VMi)g - f{u2,t2) + m(m - 1 )m™“^(Vm2, VM 2)g||xo 

< \\fiui,ti) - fiu2,t2)\\xo +m{m- 1)||m™“^(Vmi, VMi)g - M™“^(Vm 2, VM2)g||xo 

< Wfi^uti) - f{u 2 ,ti)\\xo + ||/(M 2 ,tl) -/(M2,t2)|Uo 

+m{m - l)\\uf^-^{Vui,Vui)g - ulf-^{Xui,Xui)g + u^-^{Xui,Vui)g 

(6.21) -M™-2{VM2, VMi)g + M™-2(VM2, VMi)g - M™-2(VM2, VM2)g||xo, 


with ti,t 2 G (0,To). By assumption, 

/(M,i) = ^y /(-A,t)(M +A)"^dA, 


uGU, t e [o,ro]. 
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Therefore, 

(6.22) f{ui,t) - f{u 2 ,t) = {u 2 -ui)^ J f{-X,t){ui + X)~^{u 2 + X)~^dX, fe[0,To], 
and 

(6.23) f{u,ti)-f{u,t2) = -^J {f{-X,ti) - f{-X,t2)){u +Xy^dX, uGU. 

Equation (I6.22L the embedding (I6.3p together with Lemma [6.21 Lemma [6.31 the boundedness 
of U and the fact that embeds into a Banach algebra of multipliers on Xq shows that the 

first term on the right hand side of Equation (16.211) is bounded by ci||ui — ^211x1/,, for a suitable 
constant ci. The second term is bounded by C 2 |ti —^ 2 ! in view of (I6.23P and the Lipschitz continuity 
of /(—A, t) in t. 

Neglecting the factor m{m — 1) we estimate the third term by 

(6.24) IK-2 - K-^IU(Xo) II (Vui, Vui),||xo 

+ IIK“"||£(Xo) (||(V(ui -U2),Vui),||xo + ||(V(ui -U2),Vu2)sl|xo). 

Let i be a smooth function on B that is equal to a; in [0,1 /2] x 9B and constant 1 outside [0,1] x 9B. 

Choose 0 < ^ < 1 — 2/q such that 7 + 1 — 2/q + ^ > (n + l)/2, cf. (11.411 . Given ui, U 2 € Xijq^q 
we find for the local partial derivatives that x^dvi € 'H®+i“ 2 /<j-e,(n+i )/2 e > 0 and 

x~^dv 2 G Xg, where, in the collar neighborhood of the boundary with coordinates (x,y) we use 
dx and , j = 1 , ... ,n. Hence the norm of (Vui, 'Vv 2 )g in Xg can be estimated with Corollary 
13.31 by the norms of the x^dvi as multipliers on Xg and the norms of x~^dv 2 in Xg. Thus 

(6.25) IKVui, VU2)5||xo = ||(i^Vui,X-^VU2)s||xo < Cl||ui||x,/,.Jk2||x,,,„ 
for a suitable constants cq , ci. 

This enables us to estimate the first term in (I6.24p with the help of Equation (16.191) and the 
second and third with (16.251) by caHui — U 2 ||xi/,, for a suitable constant cg and thus completes 
the argument. □ 


7. The Equation in Spaces with Asymptotics 

In this section we will obtain more precise statements on the asymptotics of the solutions to the 
porous medium equation near the tip. As before, let s G K and 1 < p,q < 00 with 

(7.1) —1 + ^ ^ + - < min{0,s}. 

P q 

For the time being, we continue to assume that 7 satisfies (imi) with 7 > ^ + I and denote by 

A, : T>(A J = ^ c ^ ny(M) 

the unbounded operator in (B) associated with these data. 

Below, we will consider the Laplacian as an unbounded operator in Yg = (D(Ag) with domain 
Yi = D^Ag). In order to obtain once more a maximally regular solution of the porous medium 
equation we will have to suppose that dim(B) ^ 3 , to make further assumptions on the spectrum 
of Ag, and to slightly change 7. We start with the following elementary result. 
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Lemma 7.1. Let 0 < a < u be a strictly positive function in 

X, ■= © C,K’'^(B)) 1 

-,g F F q^H 

We let 

(7.2) J = I1(AJ = ■H^+2’'^+2(]B) © C 
and define recursively for k G N 

:= {v e Vi{uA,)'^)\uA,v G V{{uA,)^)}. 
Then, for any 0 G [0,7r), there exists some c > 0 such that 

(7.3) c-uA,-. P((uAj'=+i) ^ V{{uA,f) 
is R-sectorial of angle 9. 


Proof. Definition CM makes sense, since m is a multiplier on by our choice of q. As 

c — uAg is sectorial on 'D{A,,), it is also sectorial in (17.111) for each k by Lemma V.1.2.3 in [5]. For 
the i?-sectoriality we proceed by induction starting with Theorem 16.11 for fc = 0. Assume that 
the assertion holds for some fc > 0. Then, for any Ai,..., Xk G {z G C \ |argz| < 9}, K G N, and 
xi,...,XK G T>((uAJ'=+i), 


K 


L2(o,l;X>((nAJfc+i)) 


^Ap(c-MA^ + Ap) 

p=i 

K 

< Ci|| '^{c - uAf){Xp{c - uA^ + Ap)“^epXp) 
p=i 
K 

^ Ap(c - uA^ + Ap)-^ep(c - uAf) 


L2(o,1;D((«AJ'“)) 


= Cl 

' P =1 
K 

- C 2 II ^ep(c-MAJ 

p=i 
K 

- y^Cpa^p 

p=i 


L2(0,1;I5((uAJ'=)) 


L2(0,l;r>((uAJ'')) 


L2(0,1;D((mAJ''+!)) 


for appropriate constants ci, C2 and C3 and c > 0. 


□ 


7.1. The domain of the bilaplacian. Associated with the extension A^ of the Laplacian given 
in (14.51) we define the bilaplacian A^ as the unbounded operator in TLp'^ (B) with domain 

V{AI) = {uGV{A,)\A,uGV{A,)}. 

According to Lemma 4.3 in 

(7.4) V{Al) = © 0 fp © c. 

p 

Here, the summation in the direct sum is over all poles of the inverse of the conormal symbol of A^ 
such that the associated asymptotics space £p is a subset of )D(Ag). In more detail; The inverse 
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of the conormal symbol has been computed in [531 Section 4.2]: 


ctm(A^)(z) 


= E' 

j=0 


{z-qj){z-q, ){z + 2-q+)iz + 2-q. ) 




where the q^ are as before and ttj is the L^-orthogonal projection onto Ej. The poles therefore lie 
in the points p = qf and p = — 2, j = 0,1,.... The asymptotics spaces have the form 


(7.5) £p = {x ^ log{x)oj{x)ei{y) + X ^uj{x)eo{y)}. 

Here, a; is a cutoff function and, for p as above, cq, ei belong to the eigenspace of the boundary 
Laplacian to the eigenvalue p; the log-terms only occur if p is a double pole. 


7.2. 

(7.6) 


Geometric Assumption. We will now suppose that 


Ai < 



n = 1 
n > 3. 


This implies that the value e„ introduced in Equation dm is larger than 3/2 for n 
than 1 for n > 2 . 

We can then restrict the choice of the weight 7 further by asking that 


(7.7) 


— - < 7 < min {£1 — 2 , 0 } for n = 1 , 


and 


(7.8) 


n — 3 

- < 7 < min 

O ' 


n — 3 
2 


+ En ~ 1, 



for n > 3. 


1 and larger 


This has an important consequence on the asymptotics spaces arising in the domain of the 
bilaplacian. Recall that a function x~p log^ x e{y) for 0 7 ^ e € (^““(^IB) belongs to 'Hp’‘^(]B) (for any 
choice of s G R and 1 < p < 00 ) if and only if Rep < {n + l)/2 — tr. In order to have £p C T’(Ag) 
we need to have £p C and thus 


Rep < 



n — 3 
2 


- 7- 


On the other hand we know that the only values of p that can arise are those of the form q^ and 
q^ — 2. This leads to the following statement. 


Lemma 7.2. Let n^2. With the choice of the weight 7 in (17.711 / (17.81) we have 

^^(A^) C ^ "^(B) 0 C for every e > 0. 

Note that + e„ > 7 + 3 6?/ the choice of 7 in (17.71) / (17.71) . 

Proof. In view of the specific form (1701) of the domain, we have to show that there is no p with 
Rep in the interval (—En, — 7 ]- 

Let us consider first the case n = I, where ei > 3/2. By (17.7L — 7 < —1/2. On the other 

hand, we have q^ = 0 , thus q^ — 2 = — 2 , q^ = ±ei, 

hence q^ — 2 > —1/2. So none of the p lies in the critical interval. 

For n > 3 we have e„ > I and — 7 < 0. 

As for the possible asymptotics, we have g/' = 0 , sog /"—2 = — 2 ,g)/ = n — 1 and thus 
gd — 2 = n — 3 > 0, gj" = —e„ < —1. Again, none of the p lies in the critical interval. □ 
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Lemma 7.3. Let n ^2 and the weight 7 be chosen as in (17.711 / (17.81) . Then 

-A ©C, 

q 

for a = + £■„ — + En — 7 ) ^ and any e > 0 . 

7 > by (17.81) . we see that a > +e„ 

Proof. According to Lemma 17.21 

(p(A 2 ),p(Aj)i,, -A ©c,h;+'’^+ 2 (b) 

Interpolation for direct sums, see e.g. Proposition 1.2.3.3 in [2], yields the desired embedding. □ 

Proposition 7.4. Let n ^ 2, and 7 satisfy (17.71) / (17.81) . For every u £ V{Ag) which is strictly 
positive on B we then have 'D{{uAg)^) = X>(A^). 

Proof. In view of (17.21) . we have P((MAg)^) = {n £ 'D{A^)\uAgV £ P(Ag)}. Since u is strictly 
positive, Lemma [6l2] shows that also u~^ £ © C. As this is a Banach algebra, we 

conclude that uA^v £ 'D{A^) if and only if A^v £ V{A^). This shows the assertion. □ 

7.3. Analysis of the Porous Medium Equation. For 7 satisfying (17.71) / (17.81) let 

(7.9) Fo = ©C, 

equipped with the norm 

ll'Wl + cjlvb = ||mi ||.^a + 2 ,T. + 2 ^gj + |c|, Ml £ c £ C. 

Lemma l 3 . 2 l imi 3 lies that Yq is a Banach algebra (up to the choice of an equivalent norm). Moreover, 
Yq consists of bounded continuous functions, i.e. Yq ^ (/(B). 

In the sequel we shall use the closed extension A of the Laplacian in Yq with domain 

(7.10) V{A) = Yi := ViAl) = h;+^’'>'+^(B) © 0 fp © C 

P 

as defined in dm). From the standard properties of operators in powers scales, see e.g. Lemma 
V.1.2.3 in [2], and Theorem 14.21 we obtain the following. 

Proposition 7.5. For c ^ M<o, 0 £ [0, tt) and (j) > 0, we have c — A £ V[9) (7 

Proposition 7.6. Let n ^ 2, and 7 satisfy (17.71) / (17.8F For each 0 £ [ 0 , 7 r) and each u £ Y^jq^q 
which is striclty positive on M, c — uA : Fi —>■ Fq is R-sectorial of angle 0 for sufficiently large 
c > 0. In particular, uA has maximal L'^-regularity. 

Proof. This is immediate from Proposition 17.41 Lemma [7.II and the fact that Y^jq^q ^ X^jq^q. □ 
The restriction on Ai now allows us to control the action of V on the interpolation spaces: 

Lemma 7.7. Let n 2, and 7 be chosen according to (17.71) / (17.81) . Assume in addition that q is so 
large that 7 < ei — 2 — 1/(17 ~ 1) for n = 1 and 7 < (n — 3)/2 + e„ — 1 — l/(g — 1) for n > 3. Then, 
in local coordinates near the boundary, the operators dx and ^dy., i £ {l,...,n}, induce bounded 
maps from Yi := (Fi,Fo)i „ to Yq. 

qi^ q'H 
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Proof. According to Lemma 1731 maps Yi/q^q to TLp^^ " ^ "^(B) continuously for every 

e: > 0, with a defined there. Our assumption (HID implies that g > 2 so that s + 3 — 2/g> s + 2. 
In order to establish the assertion it suffices to check that cr—1 >7 + 2. A short computation 
shows that our assumptions guarantee precisely this. The operator dx can be treated in the same 
way. □ 


We are now ready to state the main result of this section. For better legibility we repeat the 
assumptions. 


Theorem 7.8. Let dim(B) ^ 3 and assume that the first eigenvalue Ai of the boundary Laplacian 
satisfies condition (EH). Let s > — 1 and let the weight 7 satisfy (TTTll/dTia. Choose 1 < p,q < 00 
so large that 


1 n + 1 2 'rni j 

— 1 H- 1 — < mini 0, si and 7 < 

p q 


^ 1 - 2 -^; 


n = 1 
n > 3. 


_l_7 _ 1_ 

Denote by A : D{A) = Yi ^ Yq the closed extension of the Laplacian in Yq associated with these 
data, defined in (17.101) . 

For every strictly positive initial value uq S ^\iq,q we then find some 0 < T < Tq such that the 
porous medium equation (fTTl) . EH; considered in the space L'^{0 ,T;Yq), has a unique solution 

u e L‘}{0,T-,Yi) nW^’^0,T;Yo)). 


Proof. We only have to check that Clement and Li’s Theorem 12.91 can be applied to (IH.lll . (16.21) 
with the Banach couple Yq, Yi. Since Ti/q,q embeds to 'Hp+^’'>'+^(B) 0 C, Lemma lOl together 
with Proposition 17.61 shows that the operator has maximal L'^-regularity for any strictly 

positive uo G ^ijq.q- 

Similarly as in Equation (I6.20L property (HI) follows from the fact that Fi/q,q embeds into a 
Banach algebra of multipliers on Yg. Compared to Theorem [ESI the proof of (H2) even simplihes. 
Since Yq is a Banach algebra, Lemma 17.71 implies the estimate 

||(Vui, Vi;2)g||>0 < C'lkl|l>i/,,Jk2||Yi/,,,, Vi,V 2 G Yi/q^q. 

Estimate (I6.2ip (with Xq replaced by lb) can then be continued as follows 
< Cillui — U2IIY0 + ^2!^! — fol + C'allui — U 2 ||yi/,,,- 
This completes the argument. □ 
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